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We present a tight binding theory to analyze the motion of electrons between carbon 
nanotubes bundled into a carbon nanotube rope. The theory is developed starting from a 
description of the propagating Bloch waves on ideal tubes, and the effects of intertube motion 
are treated perturbatively in this basis. Expressions for the interwall tunneling amplitudes 
between states on neighboring tubes are derived which show the dependence on chiral angles 
and intratube crystal momenta. We find that conservation of crystal momentum along the 
tube direction suppresses interwall coherence in a carbon nanorope containing tubes with 
random chiralities. Numerical calculations are presented which indicate that electronic states 
in a rope are localized in the transverse direction with a coherence length corresponding to 
a tube diameter. 
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I. INTRODUCTION 

A carbon nanotube is a cylindrical tubule formed by wrapping a graphene sheet. Single wall carbon 
nanotubes (SWNTs) can be synthesized in structures 1 nm in diameter and microns longd. There has been 
particular interest in the electronic properties of SWNTs which are predicted to exist in both conducting 
and semiconducting formscl. Remarkably, it is possible to probe this behavior experimentally by contactiiag 
individual tubes with lithographically patterned electrodes or by tunneling spectroscopy on single tubeaa. 
However, most methods for synthesizing carbon nanotubes do not produce isolated tubes; instead the tubes 
self assemble to form a hierarchy of more complex structures. At the molecular scale tubes pack together to 
form bundles or "ropes" which can contain 10-200 tubes. X-ray diffraction reveals that a bundle contains 
tubes close packed in a triangular latticelJ, and measurements of the lattice constant and tube form factor 
led initially to the suggestion that these ropes contain primarily (10,10) nanotubes, a species predicted to be 
metallico. Subsequent work has demonstrated that the ropes likely contain a distribution of tube diameters 
and chiralitiesQ. On larger scales the ropes bend and entangle, so that the macroscopic morphology of a 
carbon nanotube sample is that of an entangled mat. Carbon nanotubes are also formed in various thick 
multiwalled species which exhibit their own unique electronic behaviorO. 

The |electronic properties of isolated SWNTs are controlled by the tube's wrapping vector, curvature and 
torsionQ. However, in ropes and in multiwalled tubes the interactions between graphene surfaces is expected to 
play a major role. This is the case even for crystalline graphite in the Bernal structure. Although an isolated 
graphene sheet is a zero gap semiconductor, the small residual interactions between neighboring graphene 
sheets with the ordered A-B stacking|-sequence lead to a small overlap of bands near the Fermi energy 
and eventually to conducting behaviorEI. Graphite is an ordered three dimensional crystal for which weak 
intersurface interactions are sufficient to establish quantum coherence for electronic states on neighboring 
sheets. Thus the electrons can delocalize both parallel to the graphene sheet and perpendicular to it. 

Reepgnizing this, several groups have attempted to estimate the energy scale for similar effects in nanotube 
ropea3. Here the situation is much more delicate, since the structure of a (10,10) nanotube does not permit 
perfect registry between neighboring tubes when they are packed into a triangular lattice. Nevertheless, it is 
possible to construct a nanotube crystal, a hypothesized ordered structure in which each (10,10) tube adopts 
the same orientation, and to study its electronic properties with conventional band theoretic methodaj. 
Theoretical studies on nanotube ropes show that intertube interactions in a nanotube crystal lead to a 
mixing of forward and backward propagating electronic states near the Fermi energy. The level repulsion 
between these branches leads to suppression or "pseudogap" in the electronic density of states, on an energy 
scale estimated to be a few tenths of an electron volt. It should be noted that these effects are qualitatively 
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different from those found in grapliite where intersurface interactions lead to band overlap and thus an 
enhancement of the Fermi level density of states. 

It has not yet been possible to extend these ideas to carbon nanotube ropes which contain a mixture of 
tubes with various diameters and chiralities. A direct calculation of the electronic structure for such a rope, 
which we define as "compositionally disordered," is quite complicated since the system has no translational 
symmetry cither along the rope axis or perpendicular to it. 

In this paper we develop a tight binding theory for the coupling between tubes. In the absence of intertube 
coupling the electronic states on an isolated tube are essentially the Bloch waves of the graphene sheet 
wrapped onto the surface of a cylinder and indexed by a two dimensional crystal momentum k. The 
essence of our theory is to develop the effects of intertube interactions t(ki,k2) perturbatively. We find 
that the effects of intertube interactions in a disordered rope are quite different from what one obtains for 
a crystalline rope. In fact, we find that compositional disorder introduces an important energy barrier to 
inter-tube hopping within a rope so that intertube coherence is strongly suppressed. We are led to conclude 
that eigenstates in a compositionally disordered rope are strongly localized on individual tubes, though they 
can extend over large distances along the tube direction. Numerical results illustrating this effect will be 
presented in this paper. We believe this physics underlies the experimental observation that charge transport 
at low temperature occurs by hopping conduction in nanotube ropes and mats. 

In section II we develop the tunneling model for describing the tight binding coupling between neighboring 
tubes in a rope. In this section we derive an analytic expression giving the tunneling amplitude t(ki,k2) 
between Bloch states on neighboring tubes indexed by momenta ki and k2. In Section III we apply the 
method to study the electronic structure of a rope crystal, and show that the model reproduces well the 
results of more complete band theoretic calculations on this ordered system. In Section IV we then extend 
the method to study the low energy electronic structure in a compositionally disordered rope and analyze 
the effects of intertube interactions perturbatively. We will also present direct numerical calculations on 
a compositionally disordered rope which probe the transverse localization of the electronic states. A brief 
discussion of the relation of these results to experimental data is given in Section V. 



II. TUNNELING MODEL 



In this section we derive an effective tight binding model which describes the coupling between the low 
energy electronic states on neighboring tubes. This coupling depends on the chirality and orientation of the 
tubes. Our starting point is a microscopic tight binding model which describes the coupling of the carbon tt 
orbitals both within a tube and between tubes, 

H^Ho+Ht. (2.1) 
Ho is a nearest neighbor tight binding model describing uncoupled tubes, 

Wo = -^ ^ t.4c„,-, (2.2) 

a <ij> 

where the index a labels the tubes and < ij > is a sum over nearest neighbor atoms on each tube. Tunneling 

between tubes is also represented by 

nT=Y.Yl taifiAiCbj + H.C.. (2.3) 

<ab> ij 

In the following we will assume that iai,bj — trai.rtj depends on the positions and relative orientations of the 
TT orbitals on the i and j atoms. 

The eigenstates of Tio are plane waves localized in an individual tube. Due to the translational symmetry 
of an individual tube the eigenstates may be indexed by a tube index a and a two dimensional momentum 
k. Of course the periodic boundary conditions imposed by wrapping the graphene sheet into a cylinder will 
give a constraint on the possible values of k. In the following, we wish to express the Hamiltonian in terms 
of this plane wave basis. We will focus on eigenstates with low energy, which have k near one of the corners 
of the graphite Brillouin zone. 



A. Plane Wave Basis 



It is useful to express the eigenstates of the individual tubes in a basis of plane wave states localized on 
either the A or B sublattice. Let us first focus on a single tube. The eigenstates may be described by 
considering a two dimensional graphene sheet with periodic boundary conditions. We will find it useful to 
consider two coordinate systems for the two dimensional graphene sheet. As shown in Fig. 1, the x and y 
axes are oriented with respect to the armchair and zig zag axes of the graphene sheet. The u and v axes, on 
the other hand are oriented with respect to the tube, with u pointing down the tube and v pointing around 
the circumference. For armchair tubes these axes coincide, and in general the angle between the axes is 
equal to the chiral angle of the tube. 

Suppressing the tube index, for the moment, we let 



(2.4) 



where 77 specifies the ^ or i? sublattice and N is the number of graphite unit cells on the tube. In this basis, 
the Hamiltonian for an isolated tube may be written 

Ho = -tTT^TkC^kCBk + H.c, (2.5) 

k 

where 
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7k = Ee''''i^. (2.6) 

Here are the three nearest neighbor vectors connecting the A and B sublattice indicated in Fig. 1. At 
low energy we may focus on the points k = aK^ + q, where a = ±1, are at the corners of the Brillouin 
zone shown in Fig. 2, and £ = —1,0,1. In the u — v system, the vectors can be written as 

aK.£ — aKo{cosu!£,siiiu!i), (2.7) 

where 

+ (2.8) 



is the angle that the Fermi vector makes with the u axis. 




FIG. 1. A graphene sheet. The x and y axes are oriented with respect to the armchair and zigzag axes, while the 
u and V axes point in the directions of along and around the tube respectively. The vectors di are the three nearest 
neighbor vectors connecting the A and B sublattices and p is the vector point to the center of a hexagon. 9 is the 
chiral angle. 



We will now focus on the point Kq. For small q, 7aKo+q — ^(v^i/2)(a9a: — 1%)- Introducing a spinor 
i^r^aii — c,,QKo+q tlic Hamiltonian may then be written, 

TYo = vipl^^{aqx(Tx + qya-y)'iJJaq, (2.9) 

where v = yfit^ajl and the r/ indices are suppressed. 

The tunneling Hamiltonian may similarly be expressed in this plane wave basis. Using (2.4) the term in 
(2.3) for the bond connecting tubes a and h is 

^ E E ^r<.r.e^(''--^-''"-'-"^4„k.C.,.k„ (2.10) 
Tart kaki, 

where ^a^h label the sublattice of the lattice site Va,h- 




FIG. 2. Brillouin zone of a graphene sheet. K^, i — —1,0, 1 label the three equivalent Fermi points and 6 is the 
chiral angle of the tube formed by wrapping the sheet in the v direction. 

The sums over lattice sites may be evaluated by introducing a Fourier transform of the tunneling matrix 
element. As detailed in Appendix A, we may write 

E E^'°"'^''"^'"""^"°"^''^^'^'''^^k„+G.k.+G.c;^k„c..k., (2.11) 

GaGb VaVb 

where 

ik„.k. = / dV,dVfci(r„rb)e-*--»+^'^--S (2.12) 

and G is a reciprocal lattice vector. 

We now specialize to eigenstates in the vicinity of the Fermi points, k = aKp + q. We may express the 
sum over the G"s as a sum over equivalent K points which are related to Kq by a reciprocal lattice vector. 
In the following we will see that this sum is dominated by the Kg, Ki and K_i, which lie in the "first star" 
in reciprocal space. Since Kg ■ t — 0, the sum becomes 

(2.13) 

with 
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T{aaVaCLa\abimb) = E e^""''"'"'^'"+''"'"^""''^"^'^'''^'"'''^ia.K„,„+q„,o,K,,,+q,. (2.14) 

We shall also find it useful to express the tunneling Hamiltonian in a basis in which the bare Hamiltonian 
describing the tubes is diagonal. This is accomplished by performing a rotation in the sublattice index space 
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to make (2.9) diagonal. Specifically, for a tube with chiral angle 6, the eigenstates will have momentum 
k = aKo + q, with {q^, Qy) = q{cos6, sin 6). Equation (2.9) is then 



Ho = vi>igaqie-"'^a+ + e'«V-)V'a,. 
Using the transformation ipaq = U{a,9)tp'^g, with 

U{a,9) = e-*l"«'-^e-'^"'^', 

the Hamiltonian becomes 

Wo = Vq{lp2gR'4''aqR - ^'IqL^agL)- 

In the {R, L) basis, the tunneling matrix has the form 

T' {aar]'aCia\ab'nb^b) = U'' {aa,9a)r,'^riaT{aaVa(la\abVb(lb)U{ab,9b)r,^ri'^, 

which may be written as 

1 

T'iaaVaCla\abv'b<lb) = J2 e^""''"^"'''""'"''''"''''''i«aK.,„+q„,a,K„^+q,M,, 



"lb' 



(2.15) 

(2.16) 
(2.17) 

(2.18) 

(2.19) 



where 



M= - 
2 



fia fib* fta fib* 

fla fib* ft-a ft-b* 
J —aaJ Ctb J —(Xa^ —(Xb 



(2.20) 



and 
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(27r - LOt). 



(2.21) 



(2.22) 



B. Tunneling Matrix Elements 



For simplicity, we suppose that the matrix elements Uj for tunneling between atoms on different tubes 
depend only on the distance between the atoms and are of the form. 



dij/ao 



(2.23) 



where dij is the distance between atoms i and j. 

It is useful to introduce two dimensional coordinates which are oriented relative to the tube's axis. Let us 
define a two dimensional vector r = (u, v), where u is the distance down the tube axis, and v is the distance 
around the tube measured from the "contact line" as shown in Fig. 3. 




FIG. 3. cross section of two parallel tubes of radius R, and with a separation b. 



Suppose the two tubes have a separation b as shown in Fig. 3. Then, the distance is given by 

d{ra,nf = K-W6)' + (i?sin^-i?sinf )2 (2.24) 
+ {b + 2R- Rcos^ - Rcos f )2. 

Since the range of the tunneling interaction qq is of order .SA, while R ~ 7Aand 6 3.4A, it is useful to 
expand (2.24) for u,v < b,R, 



c^(r.,r.) = .+ ^^ + ^. (2.25) 



26 2i? 

It follows that the tunneling matrix element has a Gaussian dependence on and r;, 



t(ra,rb) = toe"*/"''e"^['''Vo''' +^7^1. (2.26) 

We may now use (2.26) to evaluate the Fourier transform of the matrix elements. The Gaussian form allows 
this to be done simply. Using the fact that the total area of the graphene sheet is given by A''^ceii = ^t^RL, 
we find 



= ^toe-^/-e-^(l''»l^+l''^l^/^e-^(^»"-'=-)=J.„„..„. (2.27) 
Acell V 47r/t 

Wc now use (2.27) to evaluate the low energy tunneling matrix elements. Due to the exponential dependence 
on |kop the sum on Ka in (2.11) will be dominated by three terms K^i in the "first star" in which |K| = 
Kq = 47r/(3a). Specifically, estimating the parameters a = 2.5 A, 6 = 3.4 A and oq = 0.5 A, the exponent of 
the last term is approximately 87r^6ao/9a^ « 2.4. Thus the next star at V^Kq will be suppressed by a factor 
of exp(— 2(2.4)) = 0.01, justifying the first star approximation. For ka(b) = <Xa{b)^ai(j}j) +1a{b) we then have 

tk„k, = iT4<.„fc,„e-3««°('=»"-'='"')', (2.28) 

with 

= ^^e-^/^oe-^^^o^oio. (2.29) 

^cell 

Using (2.27) we then arrive at a final expression for the tunneling matrix element relating eigenstates on two 
tubes. 

1 

xe-^^"°^'=--'=^-')\„(„=a„(„K„,„(,,,,+q„,,,. (2.30) 



C. Estimate of tr from the band structure of graphite 

The tunneling model described above may be used to describe the coupling between flat graphene sheets. 
Since the transverse bandwidth of graphite is well known, this allows us to estimate the prefactor tx in the 
tunneling matrix element. The coupling between two flat graphene sheets is described by the R ^ oo limit 
of the above theory. In this case, the Gaussian dependence on kav — kbv can be written as a (kronecker) 
delta function: 



«0 -^flao(fc„„-fa»)^ 

47ri? 



havki^- (2.31) 



We thus obtain 

tk„k[, = iG^k„k[,, (2.32) 



with 



t^ = ^i„e-''/-e-^^-^c. (2.33) 

For this calculation, wc find it most convenient to use the sublattice basis for the electronic eigenstates. 
Using (2.27) the tunneling Hamiltonian for two graphene sheets is then, 

'M-T= ^("^aq|a?y6q)V'L,,„qV'f>a»)6q) (2-34) 

where 

1 

T{aija<l\avb<i) ^tGj2 e^"K'-(^''+(''"-'"''^', (2.35) 
e=-i 

and Ap — Pa — Pb- For AB stacking of graphite, Ap — r, so, using expiaK • r = the only nonzero 
term is 

r(alq|a - Iq) = 3iG- (2.36) 

Thus tunneling only connects the A sublattice on the "A" sheet to the B sublattice on the "B" sheet. This 
is to be expected, since an atom on the B sublattice of the "A" sheet sits above a hexagon on the "B" sheet. 
The Hamiltonian for an AB stacked crystal of graphene planes then has the form, 

n^Y^ wV'Lqla&Cr^ + %Cry)V'saq + ^ta ^ V'LaAqV'2s+laBq + V'LaAqV'2s-laBq, (2.37) 
s s 

where s indexes the graphene sheets. This can be simplified by introducing a transformation which inter- 
changes the A and B sublattices of the graphene lattice when a — —1, followed by a transformation which 
interchanges the A and B sublattice on the odd (2s ± 1) graphene layers. The Hamiltonian then has the 
simpler form, 

7^ = ^ w^lq • aV'. + Y 1^0^-1(1 + afT,)V's'. (2.38) 
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<ss'> 



This leads to an energy dispersion 



E{q, q,) = 3tG cos bq, ± x/v^\q\^ + (Stc cosbq,)^. (2.39) 

The bandwidth for transvecse motion is then W = 12tQ. Experimentally, the bandwidth of graphite is in 
the range W = 1.2 — 1.6 eVEI. This leads to an estimate to = 0.1 eV. 

Comparing (2.29) and (2.33) we may relate the tube tunneling matrix element to that of graphite, 



For a tube with radius R = iK we find 

tr = 7.5meV. (2.41) 



III. ELECTRONIC STRUCTURE OF A ROPE CRYSTAL 



We now apply the tunneling model described above to the problem of the electronic structure of nanotube 
ropes. We begin by considering the simpler problem of an orientationally ordered crystal of (10,10) tubes. 
We then consider a compositionally disordered rope. 
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(10,10) tubes can be arranged in a triangular lattice in which each tube has the same orientation, and the 
tubes face each other via A-A coupling, B-B coupling and hexagon-hexagon coupling. For tunneling in the 
same subband, kav = kbv Again, we find it useful to use the sublattice basis for the tube eigenstates. For 
each bond, the tunneling between the pair of tubes is described by equation (2.27) with 

1 

T{avaci\a7ji,q) =tTY. e'"K-(^''+("''+'"')^). (3.1) 

For an A-A bond it is only nonzero for rja = rib ^ I. Using a matrix notation for the j] indices: 

T^^(a77aq|a77fcq) = ^tril + a,). (3.2) 

Similarly, for a B-B bond, 

T^^(a?7,q|ar;bq) = ^ril - a,). (3.3) 

For a hexagon-hexagon bond, we have Ap — 0, so 

T^^(a77,q|a?7foq) = Stra,. (3.4) 
The Hamiltonian describing the transverse motion will then have the form, 

Ti-T = ^{li + 72) + crz(7i - 72) + 2(7:^73, (3.5) 
q 

where 

7i = e^TCOsq-ai, (3.6) 
where are the three nearest neighbor vectors in the triangular tube lattice. Then, 

E{qx, q) = 71 + 72 ± Vi^l'-c - 273)^ + (71 " 72)^- (3.7) 

From the above estimate of tx = 7.5 meV, the density of states is plotted in Fig. 4, showing a pseudogap 
feature, associated with an energy of order 12iT ~ .09 eV. The energv scale of this pseudogap agrees well 
with the results of more sophisticated electronic structure calculationaS. 




IV. COMPOSITIONAL DISORDER 



A compositionally disordered rope contains a random distribution of chiral tubes. Since tubes with different 
chiralities have different periodicities, Bloch's theorem is of little use for describing the eigenstates of the 
entire rope. Nonetheless, in the absence of coupling between the tubes, we know that the eigenstates on each 
tube are plane waves. Our approach is to describe the coupling between these plane waves perturbatively. 
We begin by considering the simpler problem of the electronic structure of two coupled nanotubes of different 
chirality. 



A. Coupling between two tubes of different chirality 



The electronic coupling between two tubes of different chirality conserves the momentum along the tube 
up to reciprocal lattice vectors in either tube. As we have argued in section II, the sum over reciprocal lattice 
vectors is dominated by the terms in which k + G arc near the first star of K points. When the coupling 
between the tubes is weak, it is useful to view this as momentum conserving coupling between states located 
near the three equivalent K points. It must be kept in mind that two states Kq + q and Ki + q in the 
vicinity of different K's arc actually the same state. 

Fig. 5 shows the Brillouin zones of two nanotubcs with chiral angles 9a and Ob oriented so that the Ua and 
Ub axes coincide. Since the tunneling Hamiltonian conserves it is convenient to view the band structure 
of the pair of tubes as a function of A;„. Consider first the band structure in the absence of coupling. The 
solid bands describe the low energy states on one tube, while the dotted ones show the states on the other 
tube. Each set of bands is replicated three times, reflecting the three equivalent K points. In Fig. 6 we 
show the band structure in the vicinity of the Fermi points with the minimum momentum mismatch, for the 
uncoupled system(Fig. 6a), and for the couple one(Fig. 6b). To lowest order, the effect of the coupling is 
only important near points of degeneracy, i.e. where we have band crossing. This occurs in two cases. The 
first, which we refer to as a "backscattering gap" occurs when the right and left moving bands on each tube 
cross. The second, which we refer to as a "tunneling gap" occurs when the left moving band on one tube 
crosses the right moving band on the other tube. 



kv 




FIG. 5. Brillouin zones of two tubes of different chiralities. The zones are rotated such that the Ua and Uh axes 
coincide. In the lower part of the figure we show the low energy band structure of the two metallic tubes. Solid bands 
belong to one tube while dotted ones belong to the other one. The bands are replicated three times, reflecting the 
three equivalent K points. 

The effect of the coupling depends on two crucial energy scales: (I) the tunneling matrix element, t, which 
we estimated in (2.27) to be less than 7.5 meV, and (2) the energy mismatch 

AE = v[aa{KaeJu - ab{Kbe,)u], (4.1) 

which determines the energy at which the right and left moving bands on tubes a and b cross. In general, 
this energy mismatch depends on the chiral angles 9a and 6b of the two tubes as well as the Fermi point 
indices i and j. In Fig. 7 we show the variation of the energy mismatch AE with the tube chirality for 
all the metallic tubes with diameters between 1.2 — 1.5 nm. Tubes which are mirror images to one another 
have the same diameter and energy difference. The offset of the u momentum is taken at a zigzag tube; 
in that case an (18,0) tube. As we see the typical energy mismatch is a few hundred meV. The fact that 
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t <^ AE simplifies the problem considerably and justifies our perturbative approach. Of course it breaks 
down in special cases when AE is zero or very small, which occurs, for instance when the two tubes are 
mirror images of each other. 





FIG. 6. (a) Band structure of the two (uncoupled) tubes in the vicinity of the Fermi point with the rninimum 
momentum mismatch Ak, which defines the important energy scale AE = vAk. (b) Band structure of the coupled 
tubes near the Fermi energy. Backscattering gaps open quadratically with the tunneling strength t, whereas tunneling 
gaps are linear in t. In general, backscattering gaps open at different energies. 
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FIG. 7. A plot showing the minimum energy mismatch AE for different metallic tubes with radii lying between 
1.2 and 1.5 nm. The energy offset is taken at a zigzag tube(18, 0) 



The tunneling Hamiltonian couples left and right movers in first order, and therefore the tunneling gap 
is linear in the tunneling strength t. The gap opens at AE/2 and since t is much smaller than AE, one 
concludes that tube-tube coupling has a small effect near the Fermi energy. 

Backscattering gaps form through second order coupling between left and right movers on the same tube, 
and hence the gap is second order in the tunneling strength; Eg ~ / AE, which is of order 1 meV. 
Furthermore, because the tunneling matrix elements are not invariant under the interchange of the two 
tubes, the backscattering gap of each tube opens at a different energy. This means that the effect on the 
density of states near the Fermi energy is weakened by the wiggling of the gaps. In what follows, we present 
quantitative arguments justifying our expectations. 

We focus first on the tunneling gaps. The movers on each tube couple in first order. In general, crossing 
occurs at three different energies, as dictated by the momentum mismatch between the Fermi points of the 
two tubes(see Fig. 5). Since we are ultimately interested in the effect of tube interactions on the states 
nearest to the Fermi level, we only consider points with the lowest lying crossing. We denote these by 
aaKai and abi^bj ■ To find the magnitude of the gaps and their offsets we need to diagonalize the first order 
perturbation matrix. We thus consider the Hamiltonian which couples the right moving states on tube a 
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with the left moving states on tube b, and we diagonalize the matrix 



v(q-Q;oKai)„ T'{aaRcia\ctbLqb) 
T'*{aaRcia\oibLcib) — y(q - abKbj)u 



(4.2) 



where q is measured from the Fermi point of tube a. For aa,cxb = +1 the magnitude of the tunnehng gap is 



l^ai . ^bj 

COS sm — - 

2 2 



(4.3) 



It is centered about an energy Ej = AE/2 above the Fermi energy. Notice that the + sign in the argument 
of the exponential is due to the fact that the tubes face each other from the outside. As we see from Fig. 7, 
the average energy separation AE 300 meV, whereas the tunnehng gap £^J < 30 meV. This means that 
in a rope with a random distribution of chirahties, the opening of such gaps will have a negligible effect near 
the Fermi level. 

Now we focus on the backscattering gaps which form near the Fermi level. Since the left and right moving 

states on the same tube do not couple in first order, we use second order degenerate perturbation theory. 
We are interested in calculating the size of the resulting gap as well as the offset of the gap. In order to 
calculate these, we need to diagonalize the matrix which arises from the second order coupling between the 
states. In general, we have tunneling between all Fermi points on each tube. Since the magnitude of the 
backscattering gaps varies quadratically with the tunneling strength and inversely with the energy difference 
AE, the most effective contributions are those with the highest tunneling strength and lowest AE. This 
argument makes us only include the set of nearest K points. Therefore, we diagonalize 



v(q - aa'Kao)u + E%ji 



Tpa 



where 

<taib> 



and A, A' = R, L. 

Diagonalizing, we get 



[T'(aaAq„|a6iq6)T'*(aaA'q„|a6iq5)-T'(aaAqa|a6i?q6)T'*(a„A'q„|a6i?q6) 



(4.4) 



(4.5) 



qa=0 



rpa 
^RR 



± V I V{CI - aa'Kao)u 



^RR 



4- I E°- 



RL\ 



(4.6) 



The backscattering gap is Eg = 2|£^^^|, and it opens around an energy offset E^ = {E'^j^ + E1^)/2. We 
thus find 



pa _ ^"-T 
' vKr 



E 

<iah> 



- 5 Kao -ffo (oa sin „ +06 sin Wi,^ . ) ^ 



aa COS LJaia ~ ^^b COS CObtt 



(4.7) 



and 



9/2 



■ ^ RaoK^ (aa sin Wa(!a sin CJ^^; )^ 



Ua COS LOala 



(If, COS LLlf,(i^ 



■ COSUbtu 



(4.8) 



Eqs(4.7) and (4.8) show that, in general, the gap offset Eo is greater than the gap Eg. In addition, one 
expects that the offsets and gaps of both tubes will generally be different. This means that the gaps wiggle 
around the Fermi energy as the chiral angle is changed, leading to the conclusion that in a rope formed of a 
random collection of chiralities, the effect on the density of states around Ep is very small. 

Let us now have a closer look at the contributions of different tunneling points. In general, only one set of 
Fermi points will dominate, unless the tubes are mirror images of each other. We now argue that in a case 
when the tubes have different chiralities, it is a certain set of Fermi points that is actually important. 
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We want to understand which tubes significantly couple to each other, and for those tubes, the Fermi 
points at which the coupling is most effective. To do this, we study the quantity 

AE vKo OLa COS ijJala, - CCfc COS UJu^ 

This quantity will be dominated by the exponential factor, and in cases where different sets of points have 
nearly equal exponential contribution, the denominator will dominate. 

For nearly armchair tubes, the maximum is at tOa(b)e{e') ~ 0, i.e., around the Kq points(Ko-Ko tunneling). 
In that case, the exponential factor is approximately 1, and the denominator takes it minimum valuc(^ 7 
meV) as the two w's are closest to zero. As the tubes shift from being armchair, the denominator increases 
as the Fermi points rotate away from the u axis, thereby making the tunneling less effective. 

For tubes which arc nearly mirror images of each other, the dominant set is also the Ko-Kq. While 
moving away from the armchair region does not significantly change the exponential contribution, it makes 
the denominator bigger, hence making the coupling less important. 

For tubes which are nearly zigzag(but are not mirror images of each other), tunneling is the least effective, 
as both the exponential argument and the denominator are big. In some of these cases, Kq-Kq tunneling 
may not be the dominant one. 

We thus conclude that tunneling is most effective between tubes which are nearly armchair, and that 
the Kq-Kq tunneling is the most important one, and hence £'^^(and similar sums) are dominated by one 
term. In other words, the sum over reciprocal lattice vectors in the tunneling Hamiltonian is dominated by 
a single term. If we ignore the other terms, then there is no reciprocal lattice vector sum, and the system 
effectively has translational invariance in the direction parallel to the tubes. This "dominant Fermi point 
approximation" simplifies our problem considerably, since it allows us to assign a conserved momentum to 
each state. This will allow us to compute the band structure for an entire rope in the following section. 

B. Compositionally disordered rope 

In this section we study the electronic structure of a nanotube rope composed of tubes with a random 
distribution of diameters and chiralities. We expect that the momentum mismatch between the Fermi points 
of neighboring tubes will suppress the tunneling and lead to localization. In real ropes, we expect 2/3 of the 
tubes to be semiconducting. As indicated in Fig. 8, this will make the localization effects even stronger. To 
emphasize our point, we consider a compositionally disordered rope with only metallic tubes. 




(a) (b) 

FIG. 8. (a) A compositionally disordered metallic rope. Different gray scales indicate different tube chiralities. (b) 
A compositionally disordered rope with 1/3 of its tubes metallic. The vacant circles denote semiconducting tubes. It 
is clear that the semiconducting tubes percolate along the rope. 

To solve the problem, we employ the "dominant Fermi point approximation" introduced in the preceding 
section. In this approximation, the momentum k = Kqu + 9 is conserved by the tunneling Hamiltonian. We 
may thus write 

i <ij> 
1 o 



(4.9) 



with 
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and 
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where T'is the tunnehng matrix given by 
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The Hamiltonian may now be diagonaUzed for each k by diagonahzing a AN x 4A'^ matrix, where N is the 
number of tubes in the rope. For each k, the rrfi^ eigenstate may be described by a "wavefunction" CS;(^)) 
which is the ampHtude for the particle to be in state a, 77 on tube i. 

A portion of the band structure of the metallic rope is shown in Fig. 9a. It is clear that there is no 
significant change in the vicinity of the Fermi energy. As we have argued before the backscattcring gaps 
wiggle around the Fermi energy, thereby negligibly changing the density of states, which is shown in Fig. 
9b. Therefore, no pseudo gap develops. 
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FIG. 9. (a) Low energy band structure of a compositionally disordered metallic rope, showing tunneling and 
backscattcring gaps. The latter wiggle around the Fermi energy, leading to a negligible effect on the density of states, 
which is shown in (b). 
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The extent of localization of the eigenstates may be quantitatively measured by computing the correlation 
function 

C{r^,E) = Yl IC,WnCv(j)l''5(|R^ - - r±)S{Em - E), (4-14) 

ijmaria' ri' 

at the Fermi energy, where is the position of tube i in the rope. As shown in Fig. 10, the correlation 
function decays exponentially with distance, C{r^,EF) oz e^^''^/^^, indicating that the eigenstates are 
localized perpendicular to the tube axes with a localization length ~ lOA. Thus the eigenstates are 
predominantly on a single tube. 
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FIG. 10. The correlation function C{r±, Ef) defined by eq(4.14), which gives a quantitative measure of the local- 
ization of the states on single tubes. 



V. CONCLUSION 



In this paper we have shown that the constraints of energy and crystal momentum conservation severely 
restrict the electronic coupling between carbon nanotubes. The electronic coupling between two nanotubes 
is only effective when the eigenstates near the Fermi energy have the same momentum, which requires 
that the graphene sheets of the two nanotubes are oriented parallel to one another. This only occurs 
when the two tubes are mirror images of one another. Thus, in contrast to a crystalline rope of armchair 
nanotubes, in which eigenstates are extended throughout the rope, we find that the electronic eigenstates of 
a compositionally disordered rope are strongly localized on individual nanotubes. 

This conclusion has important consequences for the transport properties of nanotube ropes. In particular, 
it provides a natural explanation of the nonlocal effects observed by Bockrath et alB in their multi-terminal 
conductance measurements. These effects can arise when different electrical leads make contact to different 
tubes within a rope, allowing the current in a tube to "bypass" an electrical lead which it does not contact. 

In the absence of impurities, the eigenstates will be localized on a single tube, but extend across the entire 
length of a tube. Scattering, either due to impurities or tube ends, will tend to localize the states in the tube 
direction. Paradoxically, by relaxing the constraint of momentum conservation such scattering will increase 
the coupling between tubes. Nonetheless, we are led to a picture of highly anisotropic localization. 

This picture may help to explain some apparently paradoxical transport data on nanotube mats. At low 
temperatures, nanotube mats are observed to ohey the three dimensional Mott variable range hopping law, 
R = i?oexp(To/r)i/4, with To of order 100 kN. If one uses the standard formula for isotropic variable 
range hopping and knowledge ofpthe nanotube's density of states, one extracts a localization length of order 
200A. By contrast, Fuhrer et al.tll have analyzed the scaling of the hopping conductivity with electric field 
and temperature R{E, T) = f{^E/T), and have argued that the localization length is much longer, of order 
6OOOA. Our picture of anisotropic localization offers a possible resolution to this discrepancy. In the simplest 
model of anisotropic variable range hopping, Tq depends on the geometric mean of the localization lengths, 
'^hile the scaling with electric field depends on the longest localization length, ^||. 

In this paper, we have developed a general framework for describing the electronic coupling between 
graphene based structures. This approach should prove useful for other problems, including the coupling 
between neighboring shells of multiwalled tubes as well as the coupling between crossed single walled tubes. 
Analysis of these problems will be left for future work. 



APPENDIX A: LATTICE FOURIER TRANSFORMS 



In this appendix, we work out explicitly the Fourier transform in section IIA. As shown in Fig. 1 the 
positions of the lattice may be written as = R -I- p -f 77T, which may be specified by a lattice vector R, 
and a sublattice index = ±1. The position of the center of the hexagon is given by p. Consider first a sum 
over lattice sites of the form 
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-^^/(r,)e-*- = -^^/(r)e^'^-nr=R+p+,x. (Al) 



Rt; 

The sum over R may be performed by introducing reciprocal lattice vectors G, 
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Defining the Fourier transform, 



we may then write 



/(rOe-'"^-* = E e^''-'*''/G+k. (A4) 
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